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Summary 

How dynamic load affects the pitting fatigue life of external 
spur gears was predicted by using NASA computer program 
TELSGE. TELSGE was modified to include an improved gear 
tooth stiffness model, a stiffness-dynamic load iteration 
scheme, and a pitting-fatigue-life prediction analysis for a gear 
mesh. The analysis used the NASA gear life model developed 
by Coy, methods of probability and statistics, and gear tooth 
dynamic loads to predict life. In general, gear life predictions 
based on dynamic loads differed significantly from those based 
on static loads, with the predictions being strongly influenced 
by the maximum dynamic load during contact. 

With the modified TELSGE, parametric studies were 
performed that modeled low-contact-ratio involute spur gears 
over a range of gear speeds, numbers of teeth, gear sizes, 
diametral pitches, pressure angles, and gear ratios. Dynamic 
loads and pitting fatigue lives were calculated. Gear mesh 
operating speed strongly affected predicted dynamic load and 
life. Meshes operating at a resonant speed or at one-half the 
resonant speed had significantly shorter lives. Dynamic life 
factors for gear surface pitting fatigue were developed on the 
basis of the parametric studies. The effects of number of teeth, 
gear size, diametral pitch, pressure angle, and gear ratio on 
predicted life were related to the contact ratio. In general, 
meshes with higher contact ratios had higher dynamic life 
factors than meshes with lower contact ratios. A design chart 
was developed for use in the absence of a computer and 
program TELSGE. An example illustrates the use of the design 
chart. 

Introduction 

Gears may fail from scoring, tooth fracture due to bending 
fatigue, or surface pitting fatigue. Scoring failure is usually 
lubrication related and can be prevented by proper lubrication 
and proper operating temperatures. Tooth fractures are usually 
caused by poor materials, improper design, or overloading 
and can be prevented by designing for bending stresses below 
the material’s maximum allowable stress. The American Gear 
Manufacturers Association (AGMA) has a standard practice 
for predicting gear surface pitting fatigue (ref. 1). The method 
assumes that infinite life results when the maximum surface 
contact stresses are less than the material’s endurance limit. 
Surface contact stress calculations may include a dynamic 


factor to account for gear dynamic loading. The AGMA 
recommends a dynamic factor of 1 for gear teeth of high 
accuracy but states that actual dynamic loads, computed or 
measured, can be used (ref. 1). 

Gear research authorities do not completely agree on surface 
pitting fatigue. Some state that gear materials do not have 
surface endurance limits (refs. 2 and 3), as is true for rolling- 
element bearings. In 1975, Coy developed an improved model 
for the surface fatigue life of spur and helical gears, using an 
approach similar to that for rolling-element bearings (refs. 2 
to 6). This work did not, however, include the effect of 
dynamic load. 

Early contributions to gear dynamic loading were made by 
Buckingham, Tuplin, Richardson, and Attia (refs. 7 to 10). 
More recently computer-based analytical programs have been 
developed to determine gear tooth dynamic loads (refs. 1 1 to 
16). The dynamic loads of these programs depend on such 
factors as inertia and stiffness of rotating members, tooth 
spacing and profile errors, size, and speed. The loads are 
determined by solving the equations of motion of a given gear 
mesh system. 

The objective of the present study was to combine the 
dynamic load calculation procedure of Wang and Cheng 
(ref. 14) with the NASA gear life model of Coy (refs. 2 to 
6) to determine how dynamic load affects the pitting fatigue 
life of external spur gears. NASA computer program 
TELSGE, modified to include Cornell’s gear tooth stiffness 
model (ref. 17), a stiffness-dynamic load iteration scheme, 
and a pitting fatigue life analysis, was used to predict gear 
dynamic loads and life. Parametric studies using modified 
TELSGE were performed for low-contact-ratio involute gears 
with no tooth spacing or profile errors. Gear dynamic loads 
and tooth stifftiesses were calculated as a ftmction of contact 
position and speed. On the basis of the parametric studies 
dynamic life factors for gear surface pitting fatigue were 
developed as a function of speed and contact ratio. 

Analysis 

Gear Life Model 

Current theory , life model proposed by Lundberg and 
Palmgren (refs. 18 to 20) is the commonly accepted theory 
for predicting the pitting fatigue life of rolling-element 
bearings. Because the fatigue failure mechanism is similar for 
both gears and rolling-element bearings made from high- 
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strength steel, the Lnndberg-Palmgren model for bearings has 
been adapted to predict gear life (refs. 2 to 6). Reference 6 
gives the life for a 90-percent probability of survival r; of a 
single tooth on a driver or driven gear of a mesh as 

rf = -5r "“e (1) 

where 5 is a material constant based on experimental data; 
/is the tooth face width; 2p is the curvature sum at the start 
of single-tooth contact; f is the involute surface length during 
single-tooth contact; and Q is the static tooth load, normal to 
the contact. A complete list of symbols is given in appendix A. 

The life of the complete driver gear (all teeth) Li in terms 
of driver gear rotations is 

L, = N{-^% (2) 

where N\ is the number of teeth on the driver gear and e, the 
WeibuU exponent, is a measure of scatter in fatigue life. 
Experimental research on AISI 9310 steel spur gears has 
shown gear fatigue to follow the WeibuU failure distribution 
with e » 2.5 (ref. 3). 

The life of the complete driven gear L 2 terms of driver 
gear rotations is 



where N 2 is the number of teeth on the driven gear. The mesh 
life (both driver and driven gears) in terms of driver gear 
rotations is given by 

I«= + (4) 

Expanded theory.— ‘To adapt the current gear life model for 
predictions based on gear tooth dynamic loads, the tooth was 
divided into intervals (fig. 1). The use of intervals allowed 
the current gear life model to account for load and curvature 
sums varying with contact position. The complete gear tooth 
life was determined from the interval lives and methods of 
probability and statistics. The details are as follows. 

When a pair of external spur gears is in mesh (fig. 2), the 
line tangent to the base circles of both the driver and driven 
gears is caUed the line of action. The gears begin contact when 
the outside radius of the driven gear intersects the line of 
action. As the gears rotate, the contact point occurs on the 
line of action. The contact ends when the outside radius of 
the driver gear intersects the line of action. The point at which 
the pitch circles of the driver and driven gears intersect is called 
the pitch point. The distance along the line of action from the 
pitch point to the start of contact is 


= '^rl.2 - 4,2 - '>.2 sin <f> (5) 



Figure 1.— Tooth intervals of a meshing gear pair. 
Driver gear 



Figure 2.— Basic geometry of a pair of external spur gears in mesh. 


The distance along the line of action from the pitch point to 
the end of contact is 


Z 2 = - rli - rp,i sin (6) 


2 






The contact length Z is defined as 

Z = + Z2 


(7) 


Driver gear 
rotation 


Dividing the contact length into equal-size intervals of length 
Ax gives 



( 8 ) 


where J is the total number of intervals on a tooth, and 

Xi - - z\ + {i - I ) Ajc for / - 1 to 7 + 1 (9) 

where x is the contact position along the line of action. The 
value of X is negative when contact is before the pitch point, 
zero when at the pitch point, and positive when after the pitch 
point. 

The life of each interval for a 90-percent probability of 
survival is given from equation (1) by 

r,j = for 7 = 1 to J (10) 

where B and /do not change from interval to interval. Both 
curvature sum and involute length, however, change with 
contact position. 

At the contact position the radii of curvature of the 
driver and driven gears (fig. 3) are 


= rp,i sin v> + Xi 

(11) 

f^ 2 .i = fp ,2 sin <p - Xi 

(12) 


The curvature sum at the contact position is 


1 





(13) 


For the interval the average curvature sum used in the life 
model is 


2 ^ for i=j (14) 

^ 2 

The curvature sum varied slightly with contact position for 
the example gear mesh data from table I for 100 intervals on 
each tooth (fig. 4(a)). (The contact position was made 
dimensionless by dividing by the base pitch pi^.) The plot 
shows the curvature sum to be symmetric about the pitch point 
(x = 0). This was true only because the driver and driven gears 
of the example were the same size. 



Driven gear 
rotation 


Figure 3.— Curvatures of involute teeth in contact. 


TABLE I.-BASELINE DATA FOR BOTH DRIVER AND DRIVEN GEAR 


Number of teeth 36 

Diametral pitch 8 

Outside radius, cm (in.) 6.033 (2.375) 

Base pitch, cm (in.) 0.937 (0.369) 

Face width, cm (in.) 0.635 (0.250) 

Pressure angle, deg 20 

Root radius, cm (in.) 5.318 (2.094) 

Fillet radius, cm (in.) 0.102 (0.040) 

Chordal tooth thickness, cm (in.) 0.485 (0.191) 

Normal load, N (lb) 1718 (386) 

Speed, rpm 4000 

Material Steel 


The involute surface lengths of the driver and driven gears 
for the 7 * interval (for small Ax) are 


f Ax 

fi , = ( — )jc, -h Ajc tan ^ for i =j (15) 






_ / Ax\ 

VbJ ' 


Ax tamp for i = j (16) 


The involute length is a linear function of contact position (fig. 
4(b)). Equations (15) and (16) imply that rotating a gear mesh 
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1 25 50 75 100 

Interval index, j 


(a) Average curvature sum of interval. 

(b) Length of involute surface of interval. 

(c) Average normal load on interval. 

(d) Life of interval for driver gear. Tooth life, 13 200 Mrev; gear life, 3100 
Mrev; mesh life, 2400 Mrev. 

Figure 4. — Effect of contact position on gear life parameters for gear data 
from table I. 

through equal angles produces unequal involute lengths and 
thus different-size tooth intervals (as shown in fig. 1). 

By using intervals the life model considers load that can vary 
with contact position. For the interval the average load 
used in the life model is 

Qj = for i=j (17) 


The static load variation with contact position depends on 
the number of teeth in contact (fig. 4(c)). As a pair of teeth 
begin contact, the preceding pair of mating teeth are also in 
contact. This double-tooth-pair contact occurs for intervals 1 
to 41 , and it is assumed that half the applied load is transferred 
per contact. Near the pitch point single-tooth-pair contact 
occurs (intervals 42 to 59), and all the load is transferred by 
it. Toward the end of contact, double-tooth-pair contact again 
occurs (intervals 60 to 1(X)) as the following pair of mating 
teeth begin contact. As before, it is assumed that half the load 
is transferred per contact. 

The life of a complete gear tooth ry, is determined from the 
interval lives and methods of probability and statistics where 

J \-\/e 

E vr) ( 18 ) 

j=\ / 

The complete tooth life was always shorter than the lives of 
the shortest-lived intervals (fig. 4(d)). Also, intervals with 
larger applied loads had much more influence on gear tooth 
life than intervals with smaller loads. 

The tooth lives for a driver and driven gear in mesh are 
determined by the expanded life theory and equation (18). They 
are equal if the driver and driven gears are the same size. They 
are slightly different if the driver and driven gears are different 
sizes because of curvature sums and involute lengths. The 
complete gear lives and mesh life are determined, as before, 
by using equations (2) to (4) and substituting 17, for 17. 

The total number of tooth intervals was varied from 30 to 
over 4(X) to check convergence on life. Static loads were used. 
All cases predicted the same tooth life. Gear size, diametral 
pitch, pressure angle, and gear ratio were also varied to 
compare mesh lives predicted by the current and expanded 
theories. Static loads were used. The expanded theory 
predicted mesh lives a little longer than, but within 10 percent 
of, those predicted by the current theoiy for meshes with equal- 
size gears. This difference was caused by the expanded 
theory’s curvature sum variation with contact position. Thus 
for meshes with equal-size gears the curvature sum variation 
had a small effect on life. For meshes with unequal-size gears, 
however, there were greater differences in the mesh lives 
predicted by the two theories. 

Gear Tooth Dynamic Loads 

Gear tooth dynamic load model.— The contact load of 
meshing gear teeth varies as the contact point moves along 
the line of action. This is known as dynamic load. It is mainly 
caused by single- and double-tooth-pair contact transitions, 
tooth stiffness variation along the contact, and tooth profile 
deviations from true involutes (tooth profile errors). NASA 
computer program TELSGE (refs. 14 to 16) was used to 
determine gear tooth dynamic loads. The program models 
meshing gears as a pair of rigid disks connected by a spring 
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Driver gear 



Figure 5.— Dynamic model of meshing gears. 

(fig. 5). The spring stiffness corresponds to gear teeth 
stiffnesses. 

The dynamic load model uses the equations of motion 
governing the angular displacements of the driver and driven 
gears. By converting the angular movements of the disks to 
linear displacements along the line of action, and by algebraic 
manipulation, the equations of motion are represented by a 
single differential equation, where 

MeqX-h (19) 


Ceq = (23) 

where f is the damping ratio, is the equivalent stiffness 
per unit face width (discussed in die section Equivalent gear 
tooth stiffness), and is the static load per unit face width. 

The relative displacement is determined as a function of 
contact position by using a Runge-Kutta numerical method and 
solving equation (19). The dynamic load on a gear tooth is 
determined as a function of contact position by 

Pd = KX (24) 

where P^ is the dynamic load per unit face width and K is the 
combined stiffriess per unit face width (discussed in the section 
Equivalent gear tooth stiffness). Note that when X is negative, 
the teeth separate and the dynamic load is zero. Although tooth 
profile errors can be accounted for in equations (19) and (24), 
they were beyond the scope of this study. 

Gear tooth stiff ness, ^Computer program TELSGE was 
modified to incorporate the gear tooth stiffness model of 
Cornell (ref. 17), regarded as the present state of the art. The 
stiffness model consists of tooth bending as a cantilever beam, 
fillet and foundation flexibilities, and local Hertzian 
compression, all as functions of contact position. In Comell’s 
model the deflections due to bending and fillet and foundation 
flexibilities are expressed as linear functions of load, but the 
deflections due to Hertzian effects are not linear with load. 
This makes the stiffaess of a gear tooth dependent on dynamic 
load, and equation (19) nonlinear. 

Equivalent gear tooth stiffness.— Tht stiffnesses of the 
driver and driven teeth of a mesh, kx and ^2 respectively, are 
found by the methods of Cornell (ref. 17). The combined 
stiffness K for a pair of teeth in contact is 


^ • 1.1 V 11 j . 1 . 1 j- 1 * For a single pair of teeth in contact (fig. 6(a)) the equivalent 

The dependent variable X, called the relative displacement, stiffness is v o v ^ 

is the compression of the spring along the line of action, 


>< 

II 

1 

(20) 


11 

(26) 

where 


For two pairs of teeth in contact (fig. 6(b)) the equivalent 
stiffness is 

■Si = and = r*. 2^2 

(21) 



(27) 

The equivalent mass per unit face width is 


where 



ilfl + A/2 

(22) 



(28) 

The equivalent damping coefficient per unit face width Ceq 
includes the effect of viscous damping, 


^(H)_ k[-W^ 

(29) 
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(a) Single-t(x>th-pair contact. 

(b) Double-tooth-pair contact. 

Figure 6.— Gear tooth stiffness models. 


The superscript (I) refers to the first pair of teeth in contact 
and (II) refers to the second pair of teeth in contact. The 
equivalent stiffness of equation (19) varies from double-tooth- 
pair contact at the start of mesh to single-tooth-pair contact 
and back to double. 

Iteration of gear tooth stiffness and dynamic load,— 
Because of the Hertzian compression, gear tooth stiffness is 
not independent of dynamic load. TELSGE was therefore 
modified to iterate for dynamic load (fig. 7). First the static 
load is defined. As in the example (fig. 4(c)) all the load is 
transferred per contact during single-tooth-pair contact, and 
half the load is transferred per contact during double-tooth- 
pair contact. Next the combined stiffness is determined along 
the contact position by using the static load in the Hertzian 
deflection computation. Then the dynamic load is determined 
along the contact position. Next combined stiffness is 
recalculated by using the calculated dynamic load in the 
Hertzian deflection computation. Then dynamic load is 
recalculated by using the latest stiffness values. The stiffness 
and load calculations continue until the change in dynamic load 
with each iteration becomes smaller than a preset amount. 

With modified TELSGE and the example data (table I) the 
dynamic load required only four iterations to converge to 
within 0.1 percent (fig. 8). So few iterations were required 
since the Hertzian deflection was usually only 10 to 20 percent 
of the total gear tooth deflection. The variation in equivalent 
stiffness due to double- and single-tooth-pair contact transitions 
is a major excitation in the dynamic load model (fig. 9). The 



Figure 7.— Howchart of gear tooth combined stiffness-dynamic load inter- 
action scheme in computer program TELSGE. 



Figure 8.— Effect of gear tooth combined stiffness-dynamic load interation 
scheme on dynamic load for gear data from table I. 
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Figure 9.— Effect of contact position on gear tooth combined stiffness for 
gear data from table I. 


dynamic load varied appreciably from the static when the 
operating conditions of the example were used (fig. 10). The 
maximum dynamic load during contact was about 30 percent 
greater than the static load. 
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Figure 10.— Effect of contact position on gear tooth dynamic load for gear 
data from table I. 

Gear life using dynamic loads.— The expanded gear life 
model, which accounts for variations of load and curvature 
sums with respect to contact position, was incorporated in 
modified TELSGE. The dynamic loads were used in the life 
model, where 

Qi = Pd,if for 1 = 1 to 101 (30) 


(TELSGE divides the contact length into 100 intervals.) For 
the data from table I the mesh life based on dynamic loads 
was then 50 percent shorter than that based on static loads. 
The cause was the increase in maximum load during contact 
when dynamic loads were considered (fig. 10). 


Results and Discussion 

NASA computer program TELSGE, modified to include an 
improved gear tooth stiffness model, a tooth stiffness-dynamic 
load iteration scheme, and a pitting fatigue life prediction 
method, was used to perform parametric studies. Dynamic 
loads and gear mesh life predictions were performed over a 
range of gear speeds, numbers of teeth, gear sizes, diametral 
pitches, pressure angles, and gear ratios. 

Effect of Speed on Dynamic Load and Life 

Modified TELSGE was run using the mesh data in table I 
for speeds ranging from 600 to 12 (XX) rpm. At very low 
speeds the dynamic load as a function of contact position 
(fig. 1 1) resembled the static load. However, spikes occurred 
at double- to single-tooth-pair contact transitions, and at single 
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to double. As the speed increased, the dynamic load as a 
function of contact position differed appreciably from the 
static. 

The dynamic load reached a maximum at a resonant speed 
o)„ of about 8500 rpm. At speeds below resonance the 
excitation frequency from the change in equivalent stiffness 
was lower than the resonant frequency, and the dynamic load 
was basically an oscillatory load superimposed on the static 
load. This produced peak dynamic loads greater than the static 
load. At speeds above resonance the dynamic load had a 
smoother response, with peaks lower than the static. This was 
caused by the greater inertia forces at higher speeds. The 
resonant speed can be approximated by 


0 ) 


n 


cos V? 
N 



(31) 


Here although the mean equivalent stiffness ^eq varies with 
load and speed due to Hertzian effects, its influence on is 
not significant. 

For the data in figure 1 1 the maximum dynamic load during 
contact was greatest at the resonant speed (fig. 12). It was also 
greater than the static load at speeds below resonance, with 
a secondary peak at about o)/o)n = 0.5. At speeds above 
resonance the maximum dynamic load during contact 
decreased and was less than the static load above o)/o)„ » 1 .2. 

The gear mesh life as a function of speed for the mesh data 
in table I is shown in figure 13. The dynamic life factor is 
defined as 


C, = Y (32) 

where is the gear mesh life based on the expanded life 
theory and dynamic loads and is the gear mesh life based 
on the expanded life theory and static loads (as illustrated in 
fig. 4). Comparing figures 12 and 13 shows that the gear mesh 
life decreased when the maximum dynamic load during contact 
increased. This was true even though the analysis considered 



Figure 12.— Effect of speed on maximum gear tooth dynamic load during 
contact for gear data from table I. 



the load along the complete contact length. The mesh life as 
a function of speed was lowest at resonance. 

Effect of Mass, Stiffness, and Damping on Gear Life 

The mass, stiffness, and damping of a gear mesh system 
significantly affected dynamic load and life. Modified 
TELSGE was run using the mesh data in table I while varying 
the equivalent mass and keeping all other parameters the 
same. The life-speed results were identical when plotted on 
dimensionless coordinates (as in fig. 13). Modified TELSGE 
was also run while varying the equivalent stiffness ATgq and 
keeping all other parameters the same. Again, the life-speed 
results were identical when plotted on dimensionless 
coordinates (as in fig. 13). Thus the value of the equivalent 
mass or the equivalent stiffness had no effect on the life-speed 
results when plotted on dimensionless coordinates. However, 
as expected from equation (31), different values of the 
equivalent mass or the equivalent stiffness produced different 
values for the resonant speed. The equivalent mass and 
equivalent stiffness must accurately portray the gear mesh 
being modeled for the calculated resonant speed to be accurate. 

The damping force in the dynamic load model depends on 
the gear system’s viscous friction and is usually an unknown. 
Damping ratios f between 0. 1 (in eq. (23)) and 0.2 were used 
in reference 11 to correlate analytical and experimental 
dynamic load gear tests. Here damping ratios of 0.10, 0.17, 
and 0.25 were used (fig. 14). Decreasing the damping ratio 
increased the dynamic load and thus shortened the mesh life 
at speeds near the resonant speed and one-half the resonant 
speed (o)/o)n = 1.0 and 0.5, respectively). A damping ratio of 
0.17 was used in the original version of TELSGE and was 
used in this study for all other figures. 

Effect of Speed and Contact Ratio on Gear Life 

Modified TELSGE was used to predict how speed and 
contact ratio affect dynamic load and gear life. Number of 
teeth, gear size, diametral pitch, pressure angle, and gear ratio 
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Figure 14.— Effect of damping ratio on gear mesh life for gear data from 
table I. 


were varied. The driver gear data for the different runs are 
shown in table n. The different sets had basically the same 
shape while displaced upward or downward when plotted on 
dimensionless life-speed coordinates (fig. 15). In most sets the 
mesh life was shortest at the resonant speed or one-half the 
resonant speed and was significantly shorter than the life based 
on static loads at those speeds. For all sets meshes operating 
above resonance had significantly longer life when compared 
with the static load calculations. 


TABLE n.— DRIVER GEAR DATA 


[Set 1 used for baseline; shaded area indicates parameter varied from 
baseline.] 


Set 

Number 

of 

teeth 

Pitch 

radius, 

cm 

Diametral 

pitch 

Pressure 

angle, 

deg 

Gear 

ratio 

Contact 

ratio 

1 

36 

5.715 

8 

20 

1 

1.69 

2 

20 

3.m 







1.56 

3 

28 

4.445 







1.64 

4 

44 

^ 6.985 







1.73 

5 

52 

8.255 







1.76 

6 

60 

9,525 







1.78 

7 

66 

6.985 

.12 





1.80 

8 

99 

6.985 



r 



1.85 

9 

28 

4.445 

8 

14,5 i 



1.92 

10 

20 

3.175 






1.41 

11 

28 

4.445 



15 



1.46 

12 

36 

5.715 



■ 1^ 1 


r 

1.50 

13 

36 

5.715 



20 



1.75 

14 

36 

5.715 



20 


1 

1.78 


The contact ratio c, defined as the average number of teeth 
pairs in contact, is given by 



Figure 15.— Effect of speed on gear mesh life for parametric study data from table II. 
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Figure 16.— Effect of contact ratio on gear mesh life for parametric study 
data from table II. 


For a mesh with a contact ratio of 1.6, two pairs of teeth are 
in contact 60 percent of the time and one pair is in contact 
40 percent of the time. Lx)w-contact-ratio gears have contact 
ratios between 1 and 2. In the parametric studies the contact 
ratio ranged from 1.41 to 1.92. 

For the data in figure 15 the dynamic life factor was plotted 
as a function of contact ratio in figure 16 for speeds of 
0.7 and 1.0. A sixth-order polynomial curve-fit was used to 


generate the curves. At a constant normalized speed the 
dynamic life factors were about the same for meshes with 
contact ratios between 1 .4 and 1 .6 but were significantly higher 
for meshes with higher contact ratios. 

With higher contact ratios the equivalent stiffness (fig. 9) 
had a smaller duration of single-tooth-pair contact and thus 
a smoother transition of double- to single- to double-tooth- 
pair contact. This resulted in lower dynamic load factors and 
higher dynamic life factors. For the sets studied, the resonant 
speed varied with equivalent mass, mean equivalent stiffiiess, 
pressure angle, and number of teeth. 

A general design chart for the dynamic life factor of a gear 
mesh was developed from the parametric studies (fig. 17). The 
objective was to determine the dynamic life factor as a single 
simple function of speed and contact ratio to be used when 
a computer and program TELSGE are not available. The heavy 
solid line represents the best fit of the results of the parametric 
studies. For co/oj„ < 0.5 the dynamic life factor can be read 
directly from the plot by using the scale on the left. For 
aj/o)„ > 0.5 the dynamic life factor is the product of the value 
of the curve (using the scale on the right) and the contact ratio 
to the sixth power. The light dotted lines represent the actual 
results of the parametric studies and indicate the possible error 
when using the chart. An example problem given in appendix 
B demonstrates the use of the design chart. A simplified hand 
calculation of gear tooth stiffness is also given in appendix B. 
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Summary of Results 

How dynamic load affects the pitting fatigue life of external 
spur gears was predicted by using a modified version of the 
NASA computer program TELSGE to perform parametric 
studies. TELSGE was modified to include a surface pitting 
fatigue life analysis. The parametric studies modeled low- 
contact-ratio involute gears with no tooth spacing or profile 
errors. The following results were obtained: 

1. Gear life predictions based on dynamic loads generally 
differed significantly from those based on static loads and were 
strongly influenced by the maximum dynamic load during 
contact. 

2. Gear mesh operating speeds strongly affected predicted 
dynamic loads and thus gear life. In most cases studied, meshes 
operating at a resonant speed or one-half the resonant speed 
had significantly shorter lives than the life based on static loads. 
Meshes operating above resonance had significant longer lives. 


3. In general, meshes with higher contact ratios had higher 
predicted dynamic life factors than meshes with lower contact 
ratios. 

4. Damping significantly affected predicted gear mesh life 
for meshes operating at or near a resonant speed or one-half 
the resonant speed. 

5. A solution for dynamic load converged with only a few 
iterations of gear tooth stiffness and dynamic load because the 
Hertzian deflection was relatively small in comparison with 
the total gear tooth deflection. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 3, 1986 


Appendix A 
Symbols 


B material constant (2.23 x 10* N/m''^, 35 000 Ib/in. * 
ref. 21) 

Ceq equivalent damping per unit face width, N sec/m^ 
(lb sec/in. 

Cy dynamic life factor 

c contact ratio 

d distance of inscribed parabola, m (in.) 

E modulus of elasticity. Pa (psi) 

e Weibull exponent 

/ gear tooth face width, m (in.) 

/lo tooth thickness at root radius, m (in.) 

/o beam cross-sectional moment of inertia, m^ (in.'^) 
i contact position index 

J number of intervals; or polar mass moment of inertia 
per unit face width (1/2 mrl for disk), kg m (lb sec^) 
j interval index 

K combined stiffness per unit face width, Pa (psi) 

^eq equivalent stifftiess per unit face width. Pa (psi) 

mean equivalent stiffness per unit face width. Pa (psi) 
k gear tooth stiffness per unit face width. Pa (psi) 

L life for 90-percent probability of survival, Mrev 

f involute length, m (in.) 

M effective mass per unit face width, Jfrl, kg/m 
(lb sec^/in.^) 

Afgq equivalent mass per unit face width, kg/m (lb sec^/in.^) 
m mass per unit face width, kg/m (lb sec^/in.^) 
nio module, mm/tooth 

N number of teeth 

P diametral pitch, teeth/in. 

dynamic load per unit face width, N/m (Ib/in.) 

Ps static load per unit face width, N/m (Ib/in.) 

base pitch, I'Kr^/N, m (in.) 

Q normal load, N (lb) 

Q average interval load, N (lb) 

Qj tangential load, Q cos y?, N (lb) 

R radius of curvature, m (in.) 
ri, base radius, cos m (in.) 

outside radius, m (in.) 

Tp pitch radius, Nmjl = A//2P, m (in.) 


Tr root radius, m (in.) 

s displacement, m (in.) 

T torque, N m (lb in.) 

t tooth thickness at pitch radius, m (in.) 

X relative displacement, m (in.) 

X contact position, m (in.) 

Ajc interval length, m (in.) 

Y Lewis form factor 

Z contact length, m (in.) 

Z\ contact length from pitch point to start of contact, m (in.) 

Z2 contact length from pitch point to end of contact, m (in.) 

a pressure angle at root radius, deg 
7 density, kg/m^ (Ib/in. ^) 

5 beam deflection, m (in.) 

f damping ratio 

ri life for 90-percent probability of survival, millions of 
stress cycles 

6 angular displacement, rad 
Lp curvature sum, m“^ (in. '^) 

Lp average curvature sum, m~^ (in.~^) 

(p pressure angle, deg 

0 ) speed, rpm 

o3„ resonant speed, rpm 

Subscripts: 
d dynamic life 

I gear tooth stiffness-dynamic load iteration index 

i contact position index 

j interval index 

m mesh 

max maximum during contact position 
s static life 

t tooth 

1 driver gear 

2 driven gear 

Superscripts: 

(I) first pair of teeth in contact 

(II) second pair of teeth in contact 
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Appendix B 
Example Problem 


Example problem . — Determine the dynamic life factor of 
the mesh from the data given in table HI when the driver gear 
is rotating at 5000 rpm. 

Solution.— Toe pitch radii are 


Pb = 


2xrj, 1 _ 2x(6.138 cm) 


Nx 32 teeth 

From equation (33) the contact ratio is 


= 1.205 cm (0.475 in.) 


= 6.773 cm (2.667 in.) 


Vp x = ^N\mo = ^(32 teeth) ^4.233 mm/tooth X ^ 


10 mm 


Tp 2 = ^ (100 teeth)^4.233 mm/tooth x 
= 21.165 cm (8.333 in.) 


1 cm \ 
10 mm/ 


1.858 cm ^ 

c = = 1.54 

1.205 cm 


The masses per unit face width of the driver and driven gears 
can be approximated by 


mi = yxrp i = (7833 kg/m^)x( 6.773 cm X 


1 m 
100 cm 


The base radii are 

y>=(6.773 cm) cos 25*=6.138 cm (2.417 in.) 

r *,2 = (21.165 cm) cos 25* = 19.182 cm (7.552 in.) 

From equations (5) and (6) the contact lengths from the pitch 
point to the start and end of contact are 


= 112.886 kg/m (1.637 x 10"^ lb secMn.^) 

.2 


m 2 = (7833 kg/m^x 21.165 cm x 


1 m 

100 cm^ 

= 1102.337 kg/m (1.598 X 10" * lb sec^/in.^) 


z, = V(21.590 cm)2 - (19.182 cm)^ 

- (21.165 cm) sin 25' = 0.964 cm (0.379 in.) 
Z2 = V(7.196 cm)2 - (6.138 cm)^ 

- (6.773 cm) sin 25' = 0.894 cm (0.351 in.) 

From equation (7) the contact length is 
Z = (0.964 cm) + (0.894 cm) = 1.858 cm (0.730 in.) 
The base pitch is 


The effective masses per unit face width are 

Ml = = = ^mi = i (112.886 kg/m) 

rh rlx 2 2 

= 56.443 kg/m (8.185 X 10~^ lb sec^/in.^) 

M 2 = i(1102.337 kg/m) = 551.169 kg/m 
(7.990x10“”^ lb sec^/in.^) 


TABLE m.-GEAR MESH DATA USED IN DYNAMIC LIFE 
FACTOR EXAMPLE PROBLEM 


Parameter 

Driver gear 

Driven gear 

Number of teeth 

32 

100 

Outside radius, cm (in.) 

7.196 (2.833) 

21.590 (8.500) 

Root radius, cm (in.) 

6.246 (2.459) 

20.638 (8.125) 

Lewis form factor 

0.433 

0.521 

Module, mm/tooth (Pitch, teeth/in.) 
Face width, cm (in.) 

Pressure angle, deg (rad) 

Tooth thickness at pitch radius, 
cm (in.) 

Modulus of elasticity. Pa (psi) 

4.233 (6) 
6.350 (2.500) 
25 (0.436) 
0.665 (0.262) 

2.068 X 10^* (30 X 10^) 

Density, kg/m^ (Ib/in.^) 


7833 (0.283) 


From equation (22) the equivalent mass per unit face width is 

(56.443 kg/m)(55 1 . 169 kg/m) 

” (56.443 kg/m) + (551.169 kg/m) 

= 51.200 kg/m (7.424 x 10”^ lb sec^/in.^) 

Determining teeth stiffnesses by the methods of Cornell 
(ref. 17) requires the use of a computer. For this example the 
stiffness calculations will be simplified by modeling the gear 
teeth as cantilever beams of uniform strength (beams in which 
the section modulus varies along the beam in the same 
proportion as the bending moment). The pressure angles at 
the root radii of the driver and driven gears (fig. 18) are 
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Figure 18 .— Gear tooth model for simplified stiffness calculations. 


%2 = 2(20.638 cm) 


0.665 cm 
2(21.165 cm) 


-(tan 0.378-0.378) 


+ (tan 0.436 - 0.436) 

= 1.095 cm (0.431 in.) 


From reference 23 and figure 18 the distance of the inscribed 
parabola is 


d = 


hi 

6m,Y 


where Y is the Lewis form factor. Thus the distances of the 
inscribed parabolas are 


di = 




0.1 


(0.960 cm)^ 


6m^Yi 


6( 4.233 mm/tooth x 


1 cm 
10 mm 


(0.433) 


= 0.838 cm (0.330 in.) 


— 1 I ^ i^*i I ~i 

Oil = cos I — I = cos 

Vr.l/ 

= 10.670^^ = 0.186 rad 


6.138 cm 
6.246 cm 


(1.095 cmy 


= • 


1 cin 

6( 4.233 mm/toothx-^ 1(0.521) 


10 mm 


= 0.906 cm (0.357 in.) 


/ 19.182 cm\ - 

«2 - cos ^ I = 21.651 = 0.378 rad 

' 20.638 cm 


From reference 22 the tooth thickness at the pitch radius is 
related to the tooth thickness at the root radius by 


t = 2r„ 


^ + (tan a — a) — (tan <p — <p) 


Therefore the teeth thicknesses at the root radii are 


^.1 = 2 r ,, 

= 2(6.246 cm) 


(tan a\ — ai) + (tan <p — <p) 

2rp,i 


0.665 cm 
2(6.773 cm) 


-(tan 0.186 -0.186) 


+ (tan 0.436 - 0.436) 


= 0.960 cm (0.378 in.) 


The inscribed parabola in figure 18 is a cantilever beam of 
uniform strength. From reference 24 the deflection for the 
beam is 

5 = where /q = ^fhl 

jEIq 12 

The gear tooth stifftiess per unit face width is 

Q, _ 3EIq ^ Ehl 

~ 6/ 2d^f 

For the driver and driven gears, respectively, 


k\ = 


Ehl, 

%d\ 


(2.068 X 10" Pa)( 0.960 cm X 


1 m 
100 cm 




8( 0.838 cm x 


1 m 


100 cm 

= 3.886 X 10*° Pa (5.636 x 10® psi) 
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Contact position, x 

= 4.564 X 10*° Pa (6.599 X 10° psi) 

Figure 19.— Equivalent gear tooth stillness as a fimction of contact position. 


From equation (25) the combined stiffness per unit face width 
is 

(3.886 X 10*° Pa)(4.564 X 10*° Pa) 

“ (3.886 X 10*° Pa) + (4.564 x 10*° Pa) 

= 2.099 X 10*° Pa (3.040 x 10° psi) 


From equation (31) the resonant speed in terms of driver gear 
rotations is 



For this example it is assumed that the combined stiffness 
is constant with respect to contact position. During single- 
tooth-pair contact the equivalent stif&iess is K, During double- 
tooth-pair contact the equivalent stiffness is IK. The mean 
equivalent stiffness per unit face width (fig. 19) is given as 



JTeq^+ j 


ir 3-2(9) 


(2.099 X 10*° Pa) 3 


3.571 X 10*° Pa (5.172 x 10° psi) 


From figure 17 for w/oj„ = 0.70 

^ = 0.04 
c° 

and the dynamic life factor is 

Cv = 0.04 c° = 0.04 (1.54)° = 0.53 

Thus about a 50-percent decrease in life compared with that 
using static loads is predicted for this example. Note that the 
simplified stiffness model used in the example may produce 
erroneous values for the resonant speed. The mean value of 
the equivalent stiffness per unit face width for this example was 
computed, by using Comell’s method (ref. 17) and TELSGE, 
as 2.741 X 10^® Pa (3.975 X 10^ psi). This produced a 
resonant speed of 6260 rpm and a dynamic life factor of 0.56. 
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